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INTRODUCTION 


Artificial Scientific Programming (ASP) is an 
application language designed to facilitate the computer 
solution of a large class of scientif ically oriented 
problems. The salient feature of an application language 
is the minimum effort required to program a computer 
solution. 

The ASP language is implemented on an IBM 360/91 
with an MVT job scheduler. Input is accepted from cards, 
a remote input terminal system HITS, or a question -answer 
session using the 2250 display unit. The free form 
statements are translated into a Fortran IV program, 
which is optionally compiled and executed. Programs 
generated by ASP are readable and flexible enough to 
make basic alterations a simple process, because 
general numerical algorithms frequently require modification 
to solve a difficult problem. 

In this report, the ASP language elements necessary 

to solve a system of first order ordinary differential 

equations are defined. A fifth order floating point 

1 

version of Nordsieck’s method is used to obtain 
solutions of first order initial value problems. 

Nordsieck incorporates iterative self starting, and 
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automatic revision of the step length into an efficient 

general purpose method. The step length is chosen such 

that a relative measure of the truncation error per 

integration step is less than a specified bound. To 

ensure the stability of the differential equations, it 

is necessary to compute the spectral radius of the 

2 

Jacobian matrix . Consequently , the assurance of stability 
Is left as a user supplied constraint on the step length. 



ASP SYNTAX 


I? he following conventions are adhered to in the 
description of the ASP language: 

1* Statements are order Independent unless 

explicitly stated otherwise, while expressions 
within a particular statement field are 
evaluated in the order they appear. Blanks 
are ignored except for headings. 

2. Variable names , arithmetic expressions, and 
functions should be written following the 
normal Fortran conventions, with the restriction 
that names be no more than four characters . 

3. A lower case e represents a Fortran arithmetic 
expression, c a floating point constant, k an 
integer constant, and v a variable name. 

The first statement must be one of the three listed 
below: 

30LVi£j ODiS: 

PROGRAM; DECK; 0D2: 

PROGRAM; HITS; QD3: 

For each of the above a Fortran program is generated 
which will solve a system of first order ordinary 
differential equations. If the first form is used the 
resulting Fortran program will be compiled and executed. 
Two and three will have the program punched on cards, or 
written on a disk so that it may be moved into a HITS 





library. The differential equations to be solved must 
immediately follow the SOLVE or PROGSAK request and are 
expressed as either 

Y(l f l)»e is Y( l,2)-e2; Y(l,n)»e n i 
or 

F(l)=e 1 ; F( 2)=e 2 ; F(n)=e n ; 

where Y(l,i) represents the first derivative of with 
respect to T. Che generated program is written to solve 
a system of equations with the standard nomenclature 

§ = F(T,Y) 

where Y and F are vectors. If it is more convenient to 
type in the equations with names other uhan the standard, 
the following designation must appear 
VARIABLES*. V^ V 2 ; V ? : 

For example, the system represented by 

§ = G(R.P) 

may have the ASP definition 

P( 1 , l)=e x ; P(l,2)=e 2 ; — P(l,n)=e n : 
or 

G(l)=e 1 ; G(2)=e 2 ; G(n)=e n : 

VARIABLES: R; P; G: 

If 11 computational" equations facilitate the writing of 
the differential equations, they are specified as 
EQUATIONS: v^e^; v 2 =e 2 ; V n =e n : 
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Summation.® are indicated by the ASP function 

3IGMa( v=k„ ,k ; e) 

*• 2 

To illustrate the way the EQUATION statement and the 
summation function are used consider the system of 
equations 

h = y ( vvv 

This system may be written 

Y(1,1)=Y(1)/K; Y(1,2)=Y(2)/Rj Y( 1 . 3)=Y( 3)/Hi 
EQUATIONS : S=SIGMA( 1=1 , 3: Y(I)); B*SQ3T(3)t 
SIGMA may only appear in assignment statements of the 

form v=e^+e *SIGMA( )* and its use is restricted to 

the EQUATION and DYNAMIC statements. The summand in 
the SIGMA function may contain nested sums such as the 
double sum 

S2=3IGMA(I»1,10; A( I ) *3IGMA( J*1 , 5 ; 3(I,J))) 

The size of arrays the.' appear Jn the program other 
than Y and P is transmitted by 

ARRAYS: v <k ) ; v o (k )• — v(k ), 

1 C. c. ** 

where k denotes a vector of seven or less elements. 

A convenient way to initialize arrays is to use 

DATA: v , v„, — v ; k„*c , k *c , — - k n *c : 
i * n 1 i 2 2 n n 

or 

TABULAR: list; format: 

The construction of these statements is best exemplified 
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by their ASP Fortran translations 

DATA T t . v 2 , — v n / k 1 «o 1 , k 2 *o 2 , — k n *c n / 

and 

READ (5»xx) list 
xx FORMAT (format) 

One procedure to initialize the matrices 



”i 2 3 " 


*2 3 4~ 

A « 

1 2 3 

B - 

345 


12 3 


456 


is to specify 

ARRAYS : A<3,3); B(3.3): 

DATA A: 3*1.0, 3*2.0, 3*3*0: 

TABULAR: ( (B( I , J) , 1=1 , 3) . J=1 . 3) ; 3F5-1: 

The data cards associated with each TABULAR statement 

follow the ASP program, and are arranged in the order 

requested. Statements used to represent parameters and 

Fortran arithmetic functions are of the form 

PARAMETERS: V =c ; v =c 0 ; V =c ; 

1 1 2 2 n n 

FUNCTIONS! v 1 (v,)=e ; v 0 (v„)=e,j — v (v )=e s 

i l < c c n ~n n 

The interval of integration and the print increment are 
defined in 

INTERVAL: c ; c : PRINT: c„ (defaults .0 ) : 

1 2 3 

Backward integration is implicitly defined by c^>c . 
c^ must be positive, and normally the solution will be 
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printed for T^c. , c +c , c +kc , c . 

113 13 2 

conditions are given with 

INITIAL: Y(l)^© 1 ; Y(2)*e 2 ; — - Y(n) s & n : 

Frequently, it is necessary to have functions dependent 

on the integrated values printed out along with the 

solution. This is accomplished by including 

DYNAMIC: v *e ; v --- v se : 

1 1 2 2 n n 

Another useful calculation associated with solutions 
of differential equations is the evaluation of Integrals. 
Integration is specified in general by 
INTEGRAL ( v; e) 

The limits of integration are assumed to be those in 
the INTERVAL statement. For example to evaluate 



it is necessary to include 

IMTEGRAL( A; EXP( -T) *Y( 1) ) : 

If the differential equations contain indeterminate 
forms over the interval of integration the limiting 
values are defined as 

SINGULAR: T«oj ?(!)*&, i F(2)=e 0 ; — F(n)=e : 

1 c n 

The default value for a bound on the relative truncation 
error per step is 10 . This may be changed by stating 
PRECISION: k: 

If k£5 & single precision program is written, and double 
precision if k>5. When defining a double precision 



♦ 


program, it is necessary to use the corresponding double 

precision Fortran functions if any appear in the ASP 

expressions. If a fixed step length and no automatic 

truncation error testing or step length modification 

is desired 

STEP: c : c : 

1 2 

must be included. The step length will have the value 
c^. The constant c^ is optional and if inserted 
signifies that truncation error testing is to be resumed 
when the Independent variable of Integration is greater 

than c . A heading to be printed on each page cf output 

2 

Is denoted 

HEADING: 80 or less alphanumeric characters; 

Comments may be inserted anywhere in the program and 
are enclosed by double quotes 

" any appropriate comment '» 

To have a condensed version of the output written on a 
disk so that a listing on the remote terminal is possible 
the statement 

OUTPUT: TERMINAL: 
must be included. 


» 
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'EXAMPLES 


For the first example a restricted three body 

o 

problem" (earth-moon-spaceship) is considered 


d x 

dy 


O 

x + 2 — 

- a* (x + g) - a(x - g* ) 

dt 2 

dt 


2 

d y 

dx 



y - 2— 

- a'y - ay 

dt 2 

dt 



with 

g = 1/82.45, s' = 1 - S 

a = g / ((x - g')* + y c ), a' = g' / ((x + g) 2 + y 2 ) 
and initial values 

t=0, x 0 *1.2 # x^=0 » y Q =0 , y£=-l. 04935750983. 

The solution x(t), y(t) is a closed orbit with period 
t»6. 192169331396, After a change of variable t to rewrite 
the equations as a first order system an ASP program for 
this example is 

SOLVE; ODE; 

F(l)=u(3); F(2)=u(4)j 

F(3)=U(l)+2.0*U(4)-AP*(U(l)+g)-A*(U(l)-GP)s 

F(4)=U(2)-2.0*U(3)-U(2)*(AP+A): 

VARIABLES: T; U; F: 

EQUATIONS s A=G/DSQRT ( ( U ( 1 ) -GF ) **2+U ( 2 ) **2 ) **3 : 

AP=GP/ESQRT ( ( U ( 1 ) +G ) **2+U ( 2 ) **2 ) **3 : 
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PARAMETERS! G-l .2128562765312D-2 ; 
GP=9.878?l43723469D-l! 

INITIAL: U(l)=1.2D+0; U( 2 ) =0 .0D+0 ; U( 3)=0 .ODfO ; 
U( 4)=-l .04935750983 s 
PRECISION! 9: PRINT: 6. 192169331396D-1 i 
INTERVAL! O.OD+O; 6 . 192169331396! 


HEADING: RESTRICTED 3-BODY PROBLEM: : 


Two successive colons indicate the end of the ASP 

program. The resulting solution is listed in Appendix B. 

The next example is a solution of the regular 

4 

Coulomb function which satisfies the differential 


equation 
d w 

ip 


_ L(L+1) 

(1 ' "T*— )w = 0 


for L=0,^=l/2. An ASP program to solve this equation is 
SOLVE i ODE: 


W(l,l)=tf(2); 

tf ( 1 , 2 ) = ( -1 . 0+2 . 0 * ST A/RHO-EL* ( EL+ 1 ) /RHO* *2 ) *W( 1 ) : 


VARIABLES t REOi W; F: 

PARAMETERS: EL=0.0; STA=0.5; CL=3 . 766858 E-I : 
INITIAL: H(1)=0.0; W(2)=CL: 

SINGULAR: RHO=0.0; H(1,1)=CL; tf(l,2)=CL: 

HEADING: COULOMB FUNCTION: 

INTERVAL: 0.0; 10.0: PRINT: 1.0:: 

The solution from this program is also in Appendix B. 
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APPENDIX A 

SUMMARY OF THE LANGUAGE ELEMENTS 


ARRAYS : v 1 (k 1 )j v 2 (k 2 ); --- v n (k n ): 

DYNAMIC: v^e^j ^ 2 ~ e 2 ; v n =e n J 

EQUATIONS: v^e^ v 2 =e 2 ; — v n se n l 

FUNCTIONS: ▼ 1 (l 1 )“® 1 i v 2 ^2^ =e 2 ; V n^^” e n J 

HEADING: 80 or less alphanumeric characters : 

INITIAL: Y(l)=e 1? Y(2)=e 2 ; — Y(n)=e n : 

INTEGRAL( v; e): 

INTERVAL: c, ; c_: 

1 2 

OUTPUT: TERMINAL: 

PARAMETERS: v *c ; v =c ; — v =c : 

112 2 n n 

PRECISION: k: 

PRINT: c: 

SIGMA ( v~k lt k 2 ; e) 

STEP: c^ ; c 2 : 

TABULAR: list; format: 

DATA: v lf v 2 , — v^j k 1 *c 1# k 2 *c 2# — k n *c n : 

VARIABLES: v . ; v ; v : 

1 2 3 

n any comment '* 
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APPENDIX B 
PROGRAM OUTPUT 

ASP (ODE-A) PROGRAM OUTPUT RESTRICTED 3 BODY PROBLEM ' 
PARAMETERS 

G= 1.212856276531D-02 GP= 9 . 8?87143723469D-01 

T U(I), 1=1, NE OR U(I)/F( 1) ■ 1=1, NE 

0.0 1.2000000000000D 00 0.0 

6.1921693313960D-01 9.0356308532460D-01 -5.1736905428703D-01 

1.2384338662792D 00 3.I698908458685D-OI -3 • 9249027 8842 45D-01 

1.8576507994188D 00 -4.3088791432671D-01 -5.8767318174133D-01 

2.4763677325584D 00 -9 .9989012077892D-01 -5.5615662740507D-01 

3.096084665698OD 00 -1.2624543543930D 00 7.6927307460315D-09 
3.7153015988376D 00 -9.9989013924600D-01 5.5615665735015D-01 
4.33451853197720 00 -4.30887939415940-01 5.8767326387986D-01 
4.9537354651168D 00 3 . 1698904539412D-01 3.9249016659641D-01 

5.5729523982564D 00 9 .0356307724612D-01 5.17369017336760-01 

6.19216933139600 00 1.2000000307160D 00 -3.3585912751935D-08 

-INTEGRATION STEPS = 2776- 

-TRUNCATION ERROR / 3TEP LESS THAN 1.0D-09 

ASP (ODE-A) PROGRAM OUTPUT COULOMB FUNCTION 

PARAMETERS 

EL= 0.0 ETA= 5.000000E-01 

CL= 3.766858E-01 

RHO W(I), 1=1, NE OR W(I)/F(I), 1=1, NS 

0.0 0.0 3.766858E-01 

l.OOOOOOE 00 5.166003E-01 5.929234E-01 

2.000000E 00 1.021117E 00 3.296OI8E-OI 

3.000000E 00 1.043202E 00 -3 . I69894S-OI 

4.000000E 00 4.192396E-01 -8.667134E-01 

5.000000E 00 -4.904488E-01 -8.331279S-OI 
6.000000E 00 -1.028593E 00 -1 .643893E-01 

7.000000E 00 -7.674270E-01 6.531553E-01 

8.000000E 00 1.035077E-01 9.621581E-01 

9.000000E 00 8.880071E-01 4.885538S-01 

l.OOOOOOE 01 9.391693E-01 -3.957638E-01 

-INTEGRATION STEPS = 122- 

-TRUNCATION ERROR / STEP LESS THAN 1.0E-05 


1 
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APPENDIX C 
JCL AND DECK SETUP 


Job card 

//JOBLIB DD DSN=SY3 1 . SNOBOL > DI SP=SHR 
// EX SC ASP 
//ASP. INPUT DD * 


ASP program 


/* 

//G0.DATA5 DD * 


data cards specified by TABULAR statements 
arranged In the order requested 
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